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■ Abstract. We consider stochastic spin-flip dynamics for: (i) monotone discrete surfaces in 1? 

with planar boundary height and (ii) the one-dimensional discrete Solid-on-Solid (SOS) model 
£h . confined to a box. In both cases we show almost optimal bounds 0(L 2 polylog(L)) for the mixing 

, ^ ' time of the chain, where L is the natural size of the system. The dynamics at a macroscopic scale 

should be described by a deterministic mean curvature motion such that each point of the surface 
feels a drift which tends to minimize the local surface tension [17]. Inspired by this heuristics, our 
(N ■ approach consists in bounding the dynamics with an auxiliary one which, with very high proba- 

bility, follows quite closely the deterministic mean curvature evolution. Key technical ingredients 
• are monotonicity, coupling and an argument due to D. Wilson [18] in the framework of lozenge 

tiling Markov Chains. Our approach works equally well for both models despite the fact that 
their equilibrium maximal height fluctuations occur on very different scales (log L for monotone 
i-^h ■ surfaces and \fL for the SOS model). Finally, combining techniques from kinetically constrained 

spin systems [2] together with the above mixing time result, we prove an almost diffusive lower 
bound of order l/i 2 polylog(i) for the spectral gap of the SOS model with horizontal size L and 
unbounded heights. 
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1. Introduction 

Understanding the dynamical behavior of interfaces undergoing a stochastic microscopic evo- 
lution and the emergence of mean curvature motion on the macroscopic scale is a fundamental 
problem in non-equilibrium statistical mechanics [17, 7]. Even the simpler question of rigor- 
ously establishing the correct time-scale for the relaxation to equilibrium is in many cases a 
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challenge. Similar questions arise also in combinatorics and computer science when the inter- 
face configurations can be put in correspondence with the dimer coverings of a planar graph: 
the main focus there is to evaluate the running time of Markov Chain algorithms which sample 
uniformly among such combinatorial structures (cf. in particular [18, Section 5] for background 
and motivations in this direction). In this paper we address this question for two natural and 
widely studied models and obtain essentially optimal bounds on the equilibration time. 

The first classical example is the continuous-time single spin-flip dynamics of discrete mono- 
tone surfaces with fixed boundary, cf. for instance [11, 18]. Monotone surfaces can be visualized 
as a stack of unit cubes centered around the vertices of I?, which is decreasing in both the x 
and the y direction (see Figure 1); see Section 2 for the precise definition. This is equivalent to 



Figure 1. An example of a monotone surface with planar boundary conditions 
on the plane x + y + z = 0. The picture is taken from [8] . 

the zero-temperature dynamics of the three-dimensional Ising model with boundary conditions 
that enforce the presence of an interface; alternatively, it can be seen as a stochastic dynamics 
for lozenge tilings of a finite region of the plane, or of dimer coverings of a finite region of the 
honeycomb lattice. The invariant measure is uniform over all discrete surfaces compatible with 
the monotonicity constraints and with the boundary conditions. The conjectured behavior of 
the mixing time T m i x is of order L 2 log L if L is the linear size of the region under considera- 
tion. The monotonicity condition produces dynamical constraints which prevent the application 
of standard tools to obtain non-trivial bounds on T m i x . To overcome this difficulty, a modi- 
fied "non-local" version of the dynamics, whose moves involve adding or removing piles of unit 
cubes stacked one on top of the other, was introduced in [11]. Using this device, a polynomial 
(in L) upper bound on T mix was proven in [11]. An important breakthrough was obtained by 
D. Wilson in [18], where the mixing time of the non-local dynamics was sharply analysed and 
shown to be of order L 2 log L (from below and above) . Via classical comparison arguments, this 
implies [15] that T m ; x = 0(L 6 ) for the single-site dynamics - here and below we use the nota- 
tion 0(L P ) for any quantity that is bounded above by L p up to polylog(L) factors. An improved 
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comparison, relying also on the so-called Peres-Winkler censoring inequalities [13], shows that 
T mix = 0(L 4 ); see [5, Section 4.1]. 

The second example is the continuous-time stochastic one-dimensional SOS model, described 
by a set of integer- valued heights 771 , . . . , rjL such that each rji is confined in an interval whose 
size is of order L (see Section 3 for the precise definition and Figure 2 for an illustration) and the 
heights rjo,rjL+i are fixed boundary conditions. The dynamics involves local moves where the 



(b) 





Figure 2. (a) A SOS configuration 77 on the interval [l,L], with boundary con- 
ditions and h and heights confined between — L and L + h (proportions are 
distorted for graphical convenience), (b) A SOS configuration constrained to lie 
above the "wall" fj, defined as fji = [ih/(L + 1)J (i.e. the integer approximation 
of the straight line joining the two boundary heights) . 



discrete heights can change by ±1 at each step. The invariant measure is the Gibbs distribution 
corresponding to a potential given by the absolute value of the height gradients. Due to the non- 
strictly convex character of the interaction, even obtaining a diffusive spectral gap bound (of 
order L~ 2 ) for zero boundary conditions has been a long-standing open problem. The analysis 
of an auxiliary non-local dynamics, in the spirit of [18], plus a judicious use of the Peres-Winkler 
inequalities were recently combined to obtain a mixing time upper bound 0(L 5 / 2 ) [12], while 
again the conjectured behavior is 0(L 2 log L). 

The main contribution of the present paper is a proof that, for both problems, T m i x = 0(L 2 ); 
see Theorem 2 for monotone surfaces and Theorem 3 for SOS. In the monotone surface case, 
we must restrict our analysis to the case where the boundary condition is approximately planar 
(see e.g. Figure 1), so that on a macroscopic scale, at equilibrium, the surface is flat, cf. The- 
orem 1. At the microscopic scale this corresponds to the zero temperature limit of the 3D Ising 
model with so-called Dobrushin boundary conditions, i.e. the boundary spins take the value ±1 
according to whether they lie above or below a given plane. In any finite volume, this leads to 
the uniform distribution over all monotone surfaces compatible with the given "planar" bound- 
ary conditions. In the infinite volume limit L — > 00, such uniform distribution is related to a 
translation invariant, ergodic Gibbs measure on dimer coverings of the infinite honeycomb lattice 
[9] , described by a determinantal point processes whose kernel is known explicitly. Our method 
relies crucially on the Gaussian Free Field-like fluctuation properties of such infinite volume 
state. 
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Theorem 3, together with techniques developed in the context of kinetically constrained spin 
models [2], allows us to implement a recursive analysis whose output is an almost diffusive 
lower bound gap > 1/L 2 up to polylog(L) factors for the spectral gap of the SOS model with 
unbounded heights and zero boundary conditions r/ = t)l+i = 0; see Theorem 4. 

Our approach can be roughly described as follows. At equilibrium the interface is macroscop- 
ically flat, with maximal height fluctuations much smaller than L (logarithmic in the monotone 
surface case and of order \[~L for SOS). The main step in the proof of the mixing time upper 
bounds is to show that an initial macroscopically non-flat profile approaches the flat equilibrium 
profile within the correct time. Heuristically the interface evolves by minimizing the surface 
tension and therefore feels a drift proportional to the local mean curvature. A major difficulty 
encountered in previous approaches is the possible appearance of very large (up to order L) 
gradients of the surface height. Our method consists in introducing an auxiliary "mesoscopic" 
dynamics which approximately follows the mean curvature flow, such that large gradients are 
absent in the initial condition and are very unlikely to be created at later times. The key point 
is that, thanks to monotonicity and coupling considerations, the original dynamics converges to 
equilibrium faster than the auxiliary one. 

The study of the mesoscopic dynamics involves a local analysis of the relaxation of a meso- 
scopic portion of the interface, whose size I depends i) on the current local mean curvature 1/Rt 
of the surface at time t and ii) on the size of equilibrium height fluctuations. For the monotone 

1/2 2/3 

surface case it turns out that one has to choose I ~ R t ' and for the SOS model t ~ R t ' (in both 
cases modulo polylogarithmic factors). We refer to Section 4 for a more detailed explanation of 
the strategy which leads to the choice of the different scales in the two models. Note that, since 
the evolution tends to a profile with vanishing mean curvature, R t grows with time and becomes 
much larger than the initial value Rq ~ L when equilibrium is approached. One key input is to 
prove that the equilibration time within such mesoscopic regions is of the correct order 0(£ 2 ): 
this result can be obtained following recent important progress in [5] for monotone surfaces 
and in [12] for SOS (both use the Peres-Winkler inequalities and Wilson's argument [18] in an 
essential way). The other key point, which is the main contribution of this work, is to show 
that the mesoscopic dynamics is dominated with high probability by the deterministic mean 
curvature evolution of a macroscopic smooth profile with the appropriate boundary conditions. 

In this paper we do not focus on mixing time lower bounds. Let us however mention that, 
for the monotone surface dynamics, under natural assumptions mentioned in Remark 1 below, 
one can apply [5, Theorem 3.1] to obtain a lower bound of order L 2 / log L for the mixing time. 
Similarly, for the SOS model it is known that the spectral gap is at most of order L~ 2 [12], which 
by standard inequalities directly implies a lower bound of order I? for T m i x . 

We believe our method could potentially work for a wide range of stochastic interface dynam- 
ics models where mean curvature motion is expected to occur macroscopically. An example that 
comes naturally to mind is the dynamics of domino tilings of the plane, for which at present 
only non-optimal polynomial upper bounds on T m i x are available [11]. 

A real challenge is on the other hand to prove that T m ; x = 0(L 2 ) for the monotone surface 
dynamics when the boundary height is not approximately planar, in which case the equilibrium 
shape is not macroscopically flat and arctic circle-type phenomena can occur [6]. While in 
principle our idea of mesoscopic auxiliary dynamics could be adapted to this case too, what is 
missing here are precise, finite-L equilibrium estimates on height fluctuations and on the rate of 
convergence of the equilibrium average height to its macroscopic limit. 



5 



Concerning the SOS model, a big challenge is the analysis of the dynamics for the model in 
dimension 2 + 1, for which not even crude polynomial bounds on T m i x are available, while on 
general grounds one can expect 1 once more the "diffusive" scaling L 2 log L. 

1.1. Generalities and notation. Let us recall some standard definitions for continuous time 
reversible Markov chains (see e.g. [10]). We will mostly work in the case where the state space 
17 is finite and the Markov chain is irreducible. In particular, there is a unique reversible invariant 
measure n. For £ £ 17 and t > 0, fif denotes the law of the configuration at time t started from 
the initial configuration £. The law of the chain is denoted by P. 

Given two laws /x, v on $7, we let — v\ = sup AC Q |//(^4) — v{A)\ denote their total variation 
distance. The mixing time T m j x , defined as 

T mix = inf{t>0:sup||/4-7r|| < (2c)" 1 }, (1.1) 

measures the time it takes for the dynamics to be close in total variation to equilibrium, uni- 
formly in the initial condition. It is well known that 

sup ||/4 - vr|| < e -Li/T mix J ; q 2 ) 

ten 

i.e., the worst-case variation distance from equilibrium decays exponentially with rate l/T m j X . 

If C denotes the infinitesimal generator of the reversible Markov chain, the spectral gap is 
defined as the lowest nonzero eigenvalue of — C. Equivalently, if £(f,g) = vr[/(— Cg)\ denotes 
the associated Dirichlet form, one has 

gap = irif , d-3) 

/ Var(/) 

where Var(/) stands for the variance vr[/ 2 ] — ir[f] 2 and the infimum ranges over all functions 
/ : fi i-y 1 such that Var(/) / 0. This definition makes sense also in the case where £1 is 
countably infinite, as will be the case for the unbounded SOS model to be considered in Theorem 
4 below. 

Throughout the paper, we will adopt the following conventions: 

(i) if x, y £ M. n , then d(x, y) denotes their Euclidean distance; 

(ii) if x £ R n , then we write x^ , a = 1, . . . , n for its components; 

(iii) if U C M n , then diam(U) = max{<i(x, y),x,y £ U} denotes its diameter; 

(iv) if U C IT 1 , then dU = {x £ IT 1 \ U such that 3y £ U with d(x, y) = 1}. If on the other hand 
U is a smooth subset of M n , then dU denotes its usual boundary. 

2. Monotone surfaces with "planar" boundary conditions 

Definition 1 (Monotone surfaces). A function 4> :I? ^ (Zu{±oo}) defines a (discrete) monotone 
surface if 4> x > <f> y whenever x^ < y^ a \a = 1, 2. The collection of all monotone surfaces is denoted 
byQ. 



On Q there is a natural partial order: we say that <fi < $ if 4> x < <j)' x for every x £ 

u 



Analogously, for U C I? we write <f>u < if <f> x < <f)' x for every x £ U. 



■^This model is actually quite tricky and can hide surprises: one can for instance show [4] that, if one adds a 
hard-wall floor constraint at height zero, the boundary height being also fixed to zero, the dynamics is slowed down 
by the presence of a bottleneck which causes the relaxation time to be exponentially large in L. This phenomenon, 
related to the so-called entropic repulsion, should not be present in absence of the wall. 
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2.1. Heat bath dynamics. We define a dynamics {<^(i)}t>o on monotone surfaces with initial 
condition £ and fixed boundary conditions (b.c.) outside a finite region. Let U be a finite 
connected subset of I? (the finite region) and 77 G 17 (the boundary condition) . Without loss of 
generality we will always assume that U contains the origin. 
Given <p G Q and x G Z 2 , let 0(*>+), ^ »-) g ft be defined by 



dx,+) _ J 9y 



if x^y 



min[(fe + 1, <j>( x m^ x {2)),4>( x m,xV)-l)] if X = V 



(2.1) 



and 



v 



x,~) = I 4>y if x^y 



L >0( a: ci)+i, a; C2)),0( a; ci) ja; (2) + i)] i/ s = y- 
The dynamics is a continuous-time Markov chain on the set 

ttr,,u := {(f> e n : <f> x = r) x for x £ C/}. 

The initial condition at time zero is some given £ G fi^t/- To each x G J7 is assigned an 
i.i.d. exponential clock of rate 1. If the clock labeled x rings at time t, we replace 4>{t) with 
or [</>(£)] t 1 ' - ) with equal probabilities. It is immediate to check that such Markov chain 



is irreducible and reversible with respect to the uniform measure on fi^t/, which we denote 7r^ 
or simply it. The mixing time is then defined as in (1.1) where the supremum is taken over 
£ G £l v ,U- 



2.2. Monotonicity. A function / on $7 is said to be increasing (resp. decreasing) if /(</>) < /(</>') 
(resp. f{cj)) > f((f>')) whenever <fi < cj)' . Given two laws fi, v on 0, we write fi -<v {y dominates 
stochastically /i) if < z/(/) for every increasing function /. The heath-bath dynamics is 
monotone (or attractive) with respect to the partial ordering "<", in the following sense. If /if 

denotes the law of $%{t), the dynamics at time t started from £ and evolving with b.c. rj, one has 
the following property (cf. for instance the discussion in [5, Sec. 2.1]): 

/4 ^/4,V if £<£' and t? < r/. 

In particular, letting t — > 00, one has vr^ ^ 7r^ . It is possible to realize on the same probability 
space the trajectories of the Markov chain corresponding to distinct initial conditions £ and/or 
distinct boundary conditions rj in such a way that, with probability one, 

< 4% (t) for every t > 0, if £ < £' and 77 < ?/. 

Such a construction takes the name of global monotone coupling. Throughout the paper we 
will apply several times the above monotonicity properties: for brevity, we will simply say "by 
monotonicity..." 



2.3. Mixing time upper bound. As we mentioned in the introduction, it is expected that 
T m ix = 0(L 2 logL), where L is the diameter of the region U. The next result proves such 
conjecture, up to logarithmic corrections, under the assumption that the boundary conditions 
are "approximately planar" (cf. condition (2.4) below). Such "planar" case is rather natural 
in terms of the three-dimensional Ising model: indeed, it corresponds to the zero-temperature 
limit of a system defined in the cylinder U x Z, with Dobrushin-type boundary conditions which 
are, say, "+" above some plane and "— " below. 
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Definition 2. Given neK 3 with ||n|| = 1, we write n > if > 0, i = 1, 2, 3. We let cp n G be 
the discrete monotone surface with slope n: 

Z 2 9 x H- 4>% = max{z G Z : x (1) n (1) + x (2) n (2) + zn (3) < 0} 
and n n denotes the plane 

n n = {z G M 3 : n- z = 0}. (2.3) 

For x G Z 2 we Zet n n (x) denote the vertical coordinate of the point in the plane II n with horizontal 
coordinates (x^\x^). For h > (resp. h < 0), II" is the plane obtained translating W 1 upwards 
(resp. downwards) by \h\ along the vertical direction. 

The planarity condition on the boundary conditions is specified as follows: 

Definition 3. Let n > 0. We say that rj is a good planar boundary condition with slope n if there 
exists C > such that 

|Tfc-^|<Clog(|x| + l) (2.4) 

for every x G Z 2 . 

As we mentioned in the introduction, under such boundary conditions the surface at equilib- 
rium is essentially flat: 

Theorem 1. Let n > and let r\bea good planar boundary condition with slope n. For every e > 
there exists c > such that for every a > the following holds. Let U be a finite, connected subset 
of 7? containing the origin and let L := diam(U). Then, for any L large enough, 

itl {3yeU: |0, - $\ > a(logL) 1+e ) < (1/ 'c)e- ca ^ L ^ + \ (2.5) 

We can finally formulate our mixing time upper bound: 

Theorem 2. In the same assumption of Theorem I, for L sufficiently large one has 

T mix <L 2 (logL) 12 . (2.6) 

With some technical effort (but no need of new ideas) one can improve the exponent 12 to 6 
but we will not do so, since neither is close to the conjectured optimal value 1. 

Remark 1. Concerning lower bounds: using the same idea of the proof of the lower bound on the 
mixing time of the three-dimensional zero-temperature Ising model with "+" boundary conditions 
in [5, Theorem 3.1], it is not hard to see that T m i x > L 2 /(c log L) /or a suitable c > 0/or instance 
when U = Ul^ Z 2 , with Ul cl smooth open set ofM. 2 expanded by a factor L. 

Remark 2. Concerning the assumption n > a first obvious observation is that, using lattice 
symmetries, one could replace it with the condition nW ^ 0, i = 1, 2, 3. The main reason why we 
excluded the case in which one or two components ofn vanish is that, in these cases, the fluctuations 
of the surface in U around 4> n are deterministically upper bounded by C log L, where C depends 
on the constant C in (2.4). As a consequence Theorem 1 becomes trivial and one can appeal to 
Proposition 3 below to get immediately Theorem 2. Thus, the really interesting and non-trivial case 
is n > 0. 
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2.4. Dynamics with "floor" and "ceiling". In the course of the proof of Theorem 2 we need 
an auxiliary restricted dynamics for an interface constrained between a floor and a ceiling. Let 

U and rj be as in the previous section; fix some (f> + , (jr G £l v ,u with 4>~ < cj) + and let 

nft = {<P G Q^u : <jT < < 4> + }. (2.7) 



and redefining 



ij,U 

One can define a dynamics restricted to jy simply by choosing an initial condition £ £ fi^ 



^ +) {y) = \ _i_ i /A ^ ^+1 V ^ (2.8) 

<#'-)(V) = J _ x^y 

\ max[^ - l,(l>( x m+i, x W), 0(^(1) >a; (2) + i), */ X = V 

(compare with Eqs. (2.1), (2.2)). The dynamics is again monotone in the sense of Section 2.2, 
but this time the invariant measure tt is the uniform measure onf!^. 

3. SOLID-ON-SOLID MODEL 

We turn to the study of the mixing time and spectral gap of a one-dimensional interface of 
Solid-on-Solid (SOS) type. The generic configuration (height function) of the standard SOS 
model is ?7 = (r/i, . . . , ?7l) G Z L and its equilibrium measure tt = ttlj l corresponding to boundary 
conditions 0, h G Z is 

^L,h(v) « ex p f- ^2 1 7 ^ 1 ~ ^il ( 3 - 1 - ) 

with r/o = and t)l+i = h. There is no inverse temperature parameter /3 in (3.1) since in this 
one-dimensional model its numerical value does not affect the qualitative behavior of the system 
and there is no loss of generality in fixing its value to unity. It is well known that TtL,h describes 
the law of the unique open contour in the two-dimensional Ising model in the box {l,...,L}xZ 
with Dobrushin boundary conditions (boundary spins are "+" under the line which joins (0, 0) 
to (L+ 1, h) and "— " below it), in the limit where the couplings on vertical edges tend to infinity. 

Since the mixing time (1.1) deals with relaxation to equilibrium from an arbitrary initial 
condition it is necessary to introduce the following bounded version [12] of the SOS model, en- 
closed in a rectangular box of sides of order L. Thanks to standard equilibrium estimates, see 
also Lemma 1 below, the behavior at equilibrium of this bounded version of the model is essen- 
tially the same as the usual unbounded one defined above. We come back to the unbounded 
model in Theorem 4, which deals with the spectral gap. 

For nonnegative integers L and h < L, consider the configuration space £li,h defined by 

n L ,h = {r} = (m,---,VL), m ezn[-L,L + h]} . (3.2) 

The equilibrium measure on £lL,h 1S then given by tt = ttl^ = tt(- \ £lL,h)- 

Occasionally we will consider the SOS model with further hard wall constraints, obtained by 
conditioning ttl^ to the event ^ < rj < £ 2 , where £ l G ^L,h> i = 1,2 are two configurations 
such that < £ 2 . Here, and below, we use the notation £ < a, for the natural partial order 
in flL,h defined via < a\, for alH = 1, . . . ,L. We refer to £ 1 ,^ 2 as the floor and the ceil- 

ing, respectively, and write tt\ for the corresponding equilibrium measure. If A denotes the 
maximal configuration in 0,L,h, i-e- Aj = L + h, we sometimes consider the model with £ 2 = A 
and t, 1 = fj where fji := [ih/(L + 1)J for all i = 1, . . . ,L, i.e. the interface is above the straight 



and 
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line connecting the two boundary values, cf. Figure 2 (b) . In this case one speaks simply of an 
interface above the wall. Note that fj is the SOS equivalent of the "monotone surface <fi n with 
fixed slope", cf. Definition 2. 

In the following, whenever we do not explicitly mention floor and ceiling, it is understood 
that we are talking about the bounded model where £ x = A and £ 2 = V, where V is the minimal 
configuration in flL,h- v « = 

3.1. Dynamics. The evolution of the interface is given by the standard heat bath dynamics, i.e. 
single-site Glauber dynamics described as follows. There are independent Poisson clocks with 
mean 1 at each site i e { 1 , . . . , L} . When site i rings, the height rji is updated to the new value 
maxjjjj — 1, — L} or mm{rn + l,L + h} with probabilities pi -(rj), Pi +(77) respectively, determined 
by: 

e~ 2 1 1 

Pi = i + e -2 1 {%<«} + 2 1 {6>*>«»} + l + e -2 1 {'K>6} 

e~ 2 1 1 

Pi.+ fa) = 1 + e _2 1 {^>6> + 2 1 {fe>%>"} + i + e -2 1 {%<«> (3 - 3) 

where a := min{7/j_i, 7/1+1} and 6 := max{?7j_i, r/i+i}. With the remaining probability 1 — 
(pi-{rj) + pj i+ (?/)), 77j stays at its current value. It is not hard to check that this defines a 
continuous time Markov chain with state space FlL,h an d stationary reversible measure given by 
itL,h' m the sequel we will write 77^ (rj) for the random variable describing the state of the Markov 
chain at time t with initial state £ and fi\ for its distribution. Let T m i x = T m i x (L, h) denote the 
mixing time of this Markov chain. 

We may consider the evolution of the system under hard wall constraints £ 1 , £ 2 as above. This 
amounts to the same dynamics except that any update which would violate the constraints is 
rejected. The dynamics is then reversible w.r.t. the equilibrium measure 7r| ^ associated to the 

floor £ x and the ceiling £ 2 , see also Section 2.4 above for the analogous constrained dynamics 
in the monotone surface case. In either case, with or without hard walls, the monotonicity 
considerations recalled in Section 2.2 apply here without modifications, with the natural partial 
order on configurations introduced above. 

Our main result about the mixing time of the SOS interface is: 

Theorem 3. There exists a > such that for any L sufficiently large, uniformly inO < h < L one 
has 

T mix <L 2 (logL) Q . (3.4) 
The same bound holds for the interface constrained to stay above the wall fj. 

Remark 3. 

1. It will be clear from the proof that the above bound is satisfied as soon as a > 21. With some 
effort this power of log L can be considerably improved. However, the method we use does not 
seem to be capable of reaching the presumably optimal bound T m i x = 0(1? log L). 

2. Minor modifications of the proof show that the result of Theorem 3 continues to hold as it is for 
the dynamics constrained between a floor ^ and a ceiling £ 2 for any £ 1 ,£ 2 G ^L,h sucn that 

3. A further extension is obtained by letting the constraints —L < rji < L + hbe replaced by —M < 
Vi — M + h, where M is an independent parameter, possibly much larger than L. It is possible 
to extend the proof of Theorem 3 to get the essentially sharp bound T m i x = 0(L max(L, M)). 
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3.2. Spectral gap. The unbounded version of the SOS model is given in (3.1). We write again 
7? = TTi,h for the corresponding Gibbs measure. The dynamics is the same as above except for 
the absence of the constraints V < rj < A, i.e. when the clock labeled i rings, r/i is updated to 
the new value rji ± 1 with probability p^± given by (3.1). The infinitesimal generator is given by 

L 

WW = E {Pi,+(vWi,+f(v)+PiAvWi,-f(v)} , (3-5) 

i=l 

where Vi,±f(ri) = /(r/ 4,± ) — f(r/), and rj 1 ^ is the configuration coinciding with rj everywhere 
except that at site i the value of r\i is replaced by rji ± 1. The Dirichlet form is given by 

1 L 

£(/,/) = -5]7f[p i , + (V i , + /) 2 +p i ,_(V i ,_/) 2 ] (3.6) 

i=l 

where tt[-] denotes expectation w.r.t. tt. Note that, for each finite L, C defines a bounded self- 
adjoint operator in L 2 (Z l ,tt). The associated spectral gap is defined by (1.3), where / ranges 
over all / e L 2 (Z l ,tt) with nonzero variance. 

Theorem 4. Let a > be as in Theorem 3. For some constant c > 0, for alio < h < L, the spectral 
gap of the unbounded SOS dynamics satisfies 

gap(L) > cL~ 2 (logLy a . (3.7) 

The same bound holds for the interface constrained to stay above the wall fj. 

This estimate is optimal, modulo the logarithmic factor: an upper bound 0(1/ L 2 ) is given 
e.g. in [12]. The lower bound gap(L) > c/L 2 was proven in [14] for a modified version of the 
SOS model with weak boundary couplings; such modified model is much less sensitive to the 
boundary conditions and has a genuinely different dynamical behavior. If instead the absolute 
value interaction potential were replaced by one with strictly convex behavior at infinity, then 
the correct lower bound gap(L) > c/L 2 would follow by well-established recursive methods, see 
e.g. [3]. 

4. Strategy of the proof 

As already announced in the introduction, despite the fact that the equilibrium fluctuations of 
the interface in the two models are very different, our bound T m i x = 0(L 2 ) is proved following 
a common strategy that we sketch here. 

The crucial step is the following (see Propositions 1 and 6 below for a precise formulation in 
the case of monotone surfaces and SOS model) : 

Step 1. Starting from any configuration, after time 0(L 2 ) the distance between the interface and 
the flat profile is not larger than its typical equilibrium value xl- 

At that point one can conclude T m i x = 0(L 2 ) provided that a result of the following type is 
available (see e.g. Proposition 5 below in the case of the SOS model): 

Step 2. If the initial condition £ is at distance \Lfrom the flat profile, then \\^ — tt\\ 1/or some 
T = d(L 2 ). 

The proof of such result is model-dependent: for the SOS model it was given in [12] and for 
monotone surfaces it follows from results in [5] (see Propositions 2 and 3 below) . 

In turn, Step 1 follows if one proves that, with high probability, the interface started from the 
maximal configuration stays below a deterministic interface evolution which after time 0(L 2 ) is 



11 



at the correct distance 0{\l) from the flat profile. It turns out that it is actually sufficient to de- 
fine the deterministic interface evolution along a sequence of deterministic times t n , < n < M. 
At all times t n , the deterministic interface is the boundary of C Un where, given u > 0, C u is a 
spherical cap (if we are considering two-dimensional interfaces like monotone surfaces) or a cir- 
cular segment (in the case of one-dimensional interfaces like SOS) of height u and base of linear 
size pl roughly of order L, see Figure 3. The base of C u lies on the plane/line which contains 
the macroscopic flat profile. The evolution of C u , by a kind of "flattening process", in the time 
interval (t n _i,i„] transforms C Un _ 1 into C Un . The sequence of increasing times {t n }o< n <M and 
of decreasing heights {u n }o< n <A/ will be introduced in a moment. The "domination statement" 
then is of the following type (see Propositions 4 and 7) : 

Claim 1. For alio < n < M, with high probability the following holds. For all times in [t n , L 3 ] the 
evolution started from the maximal configuration stays below the boundary ofC Un . 

The initial height uq is taken to be proportional to L and one sets to = 0; this guarantees that 
the statement of Claim 1 holds trivially for n = 0. In order to choose u n+ \ given u n , one uses the 
following procedure. Consider the spherical cap/circular segment C Un and choose a point on its 
curved boundary (e.g. the highest one). Move inward (i.e. inside C Un ) the tangent plane/line 
at the chosen point by an amount A and call c?a the diameter of the intersection between the 
plane/line with C Un , see Figure 3. 



Figure 3. The spherical cap/circular segment C Un . The deterministic evolution 
at time t n coincides with the curved portion of the boundary of C Un . At that time, 
the true stochastic evolution (wiggled line) stays with high probability below it. 
Elementary geometry shows that A ~ d 2 A (u n /p L ). The requirement u n — u n+ \ = 

A « - d A n then leads t0 U n ~ u n+l ^ (p|/Wn) 7/(2_7) . 

Then u n — u n+ \ is chosen as the critical value A n such that the equilibrium fluctuations on 
scale g?a„ are of order A„ (apart from logarithmic corrections), i.e. A n = 0(xd An )- Also, M 
is the smallest index such that um < XL, i-e- um is of the order of the equilibrium height 
fluctuations on scale L. As for the time sequence {t n }o< n <M, one sets t n+ \ — t n to be of order 
d\ n (again neglecting logarithmic corrections) : that this is the correct choice is guaranteed by 

a careful use of Step 2, applied with L = d& n . It is not difficult to realize that iju = 0(L 2 ). 
Indeed, assume for definiteness that xl ~ L" 1 for some < 7 < 1, where if 7 = we mean that 
Xl ~ polylog(L). Then, simple geometric considerations show that 




PL 



~ L 
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Approximating the recursion for u n with a differential equation gives 
since both uo and are of order L. In particular, one has roughly M = o(l(2-2 7 )/(2- 7 ))_ Then, 

M-l M-l 4/(2-7) ^ 

*M= ^(Wl-^n)^pi 2 " 7 E Z-m^) ~ L 2 7 /(2- 7 ) E L (2-2 7 )/(2- 7 ) _ " = d ^ 
n=0 n=0 n=0 

since the last sum is of order log L. Remarkably, the order of magnitude of tjw does not depend 
on the fluctuation exponent 7, while the sequence (t n , u n ) and the value of M do. The statement 
of Claim 1 for n = M allows to conclude Step 1: the evolution started from the maximal 
configuration, at time tu = 0(L 2 ), is below the deterministic evolution, which is within distance 
Xl from the flat profile. 

Another way to understand the choice of the time-scales t n is the following. If one imagines 
that the boundary of C u evolves by "mean curvature", i.e. feeling a inward drift proportional 
to the inverse of its instantaneous radius of curvature, then the time t n+ \ — t n to transform 
C Un into C Un+1 must be 0(R n x (u n — u n+ i)), where R n is the radius of curvature of C Un . One 
can easily check that, apart from logarithmic corrections, this coincides with the requirement 

t n+1 -t„~ d\ n . 

5. Monotone surfaces: Proof of Theorem 1 

An essential tool in the proof of Theorem 1 are the translation invariant, ergodic Gibbs mea- 
sures, with given slope n, on the set of monotone surfaces [9, 16]. The trick of using the 
properties of such infinite-volume states to obtain fluctuation bounds for surfaces with fixed 
boundary conditions around a finite region was already crucial in [5]. The most relevant result 
for this purpose is the the following. 

Theorem 5. [9, 16] Given n > 0, h G Z and x G Z 2 , there exists a unique law p n ,h,x on such 
that 

(0 ^n,h,x(4>x = h) = 1; 

(ii) for every x € 1?, /J, n ,h,x(<^x) = h + IL n (x — x) (cf. Definition 2); 

(Hi) for every xx,...,Xk € Z 2 and v G I? one has that the joint law of {4> Xi — (j} Xj }i<i,j<k ! ' 5 
the same as the joint law of {4> Xi +v — <j>xj+v}i<i,j<k (translation invariance of the law of the 
height gradients); 

(iv) for every finite f/cZ 2 such that x ^ U and every monotone surface £ G such that £ s = h, 
the measure fi n .h,x conditioned to {cf> x = £ x for every x ^ U} is ir^, i.e. the uniform measure 
over all 4> GO which coincide with £ outside U (DLR property). 

Concerning height fluctuations under these Gibbs measures one can prove that for every e > 
there exists a positive constant c such that, for every a > and L large enough, 

Mn,M (3y G Z 2 : d(x,y) < L and \<f> y - fi nA x(4>y)\ > a(logL) 1+e ) < V ca , (5.1) 

(see [5, Proposition 5.7]; the proof is given there for e = 1/2 and a = 1/8 but it works identically 
for e > 0, a > 0). In other words, the surface lies on average on a plane of slope n and the height 
difference between two points x, y does not differ from the average height difference by much 
more than the logarithm of d(x,y). The estimate (5.1) is obtained via the well-known fact 
[8] that height differences between two points can be written as the number of points of a 
determinantal point process whose kernel is explicitly known. 
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The connection between the measure fj, nt h,x an d the "infinite volume states" for dimer cover- 
ings of the infinite honeycomb lattice defined in [9] is well known and is discussed for instance 
in [5, Section 5]. Let us just recall that the components = 1, 2, 3 are directly related to 

the fractions of dimers of the three possible types (horizontal or rotated by ir/3 or by (2/3)tt). 

We now turn to the proof of (2.5). By symmetry it is enough to show that 

ttI (3yeU:0 y >$ + a(logL) 1+t ) < (l/c) e - c °CK*£) 1+e . (5.2) 

By monotonicity the probability in the left-hand side of (5.2) is increased whenever rj is replaced 
by rf such that r]au < Vau- m particular, let rf be sampled from fJ. n ,h,x(-\A) where x is a given 
element of dU, h is such that /i n ,M(<Ar) = Vs + ( a /2)(log L) 1+e and A is the event A = A V) u = 
{(f) e ft : 4>du > t]qu}. From (5.1) and the fact that diam(U) = L one sees that 

\\Vn,h,x(-\ A ) ~ Mn,Mll = eX P (-^ C ( lo g L ) 1+e 



3 

This is because, if y e dU, (f> y < r/ y implies <f) y < pi n> h,x{4>v) ~ ( a /^ + L) l+e , recall 

Definition 3, so that /i n ,/i,i(^ c ) = 0(exp(— (a/3) c (log L) 1+e )). Therefore, the probability in the 
left-hand side of (5.2) is upper bounded by 

J Vn,hd d v')4 {3y€U:<t>y-4$> a(logL) 1+e ) + O [exp (-| c (log L) 1+e 
= ^n,M (31/ G U : <f> y - $ > a(logL) 1+e ) + O [exp (-| c (log 

r / n . \ 1 

= O 



exp (^--c(log L) 

where in the first step we used the DLR property and in the second one the fluctuation bound 
(5.1), together with the fact that with our choice of h one has fin,h,x(4>y) = 4>y + ( a /2 + 
o(l))(logL) 1 + e . 

6. Monotone surfaces: Proof of Theorem 2 

In this section the slope n > 0, the region U and the good planar boundary condition rj 
with slope n are fixed as in Theorems 1, 2. Denote A (resp. V) the maximal (resp. minimal) 
configuration of fi^jy with respect to the partial ordering "<" and recall that <ffi (t) denotes the 
configuration at time t started from £. 

The first key ingredient is a result saying that, after time of order L 2 (log L) 12 , the surface is at 
most at distance (logL) 3 / 2 away from the plane n n of slope n (cf. Definition 2). It is here that 
the new ideas of mimicking the evolution by mean curvature play a crucial role. 

Proposition 1. Let T = (l/2)L 2 (log L) 12 . Then, there exists c > such that 

P (max(^(T) - O > (logL) 3 / 2 ) = O (V^) 372 ) (6.1) 



and similarly 



min(^(r) - C) < -(log Lf 'I 2 = O ( e ~<^ L ^ ) . (6.2) 



The second result says that once the surface is within distance (logL) 3 / 2 from n n , it does 
not go much farther than that for a time much longer than L 2 (logL) 12 . This second step is 
much more standard and its proof combines monotonicity and reversibility together with the 
fluctuation bounds of Theorem 1. 
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Definition 4. Let <£ + (resp. ) be the maximal (resp. minimal) configuration in the set 

{cP G n„ >t7 : max \<j> x - ft\ < 2(logL) 3 / 2 }. 

x£U 

Proposition 2. Let £ G Q Vj u be such that 

max|^-C| < (logL) 3 / 2 . (6.3) 

Then, there exists c > such that 

P (at < L 10 : £ 0**) = O ( e -<=0°^) 3/2 ) (6.4) 

with $7*^} given in (2.7) and <I> ± as in Definition 4. 

Finally, the last step shows that if the surface evolves constrained between a ceiling and a 
floor which are within distance 0((log L) 3 / 2 ) from IP, then mixing occurs within a time 0(L 2 ). 

Proposition 3. [5, Theorem 4.3] Let cj^ G Cl v ,u with 0~ < <p + . For the dynamics restricted to 
Q,t TJ (cf. Section 2.4) one has 

T mix = 0(L 2 (logL) 2 # 2 ) 

where H = max xe u((l>t — <f>~). 

We can now easily put together Propositions 1 to 3 to obtain the desired upper bound (2.6) 
on the mixing time: 

Proof of Theorem 2. It is a standard fact that 

max — 7r|| < L 3 max(||/4 — vr 1 1 , ||/t^ — vr 1 1 ) (6.5) 

(see e.g. [5, Lemma 6.2] for a similar statement) so it is sufficient to prove that 

max(||^-7r||,||^-^||)<l/(2 e L 3 ) 

for t = L 2 (logL) 12 = 2T. Let us consider e.g. the case of the maximal initial condition A, the 
other case being analogous. 

Define W = {0 G Q V) u ■ max xe (/ \<p x - $£| < (logL) 3 / 2 } and, for £ G 0', 

r = inf{t > : max|^|(t) - 0"| > 2(logL) 3/2 - 1}. 

x&U 

Let A be a subset of £l v ,u- Then, using Proposition 1, 

^ T (A) = i4(&(A)\Z e «') + O (e- c(logL)3/2 ) • 
Next, from Proposition 2, one has for every 

4(^4) = W{<f>H T ) G A;r > T) + O (Vc(l°g£) 3/2 ) = P $± (^(r) ei;r>T) + ( e - C (lo g L) 3 /^ 

= P $± (^(T) G A) + O (Vc(i°g^) 3/2 ) 

where P** denotes the law of the dynamics restricted to the set fi*^. Indeed, up to the random 

time r the two dynamics P and P* can be perfectly coupled so that they coincide. In particular, 
r has the same law under P and P $± . Finally, thanks to Proposition 3, T is at least (log L) 2 times 
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the mixing time of the restricted dynamics (which is 0(L 2 (log L) 5 )). Therefore, from (1.2) and 
the fact that the invariant measure of the restricted dynamics is 7r(-|$7*^}), one has 

^(/(T) £ A) — 7r(A|fi*£)| < e-^) 2 . (6.6) 
Thanks to Theorem 1 one has ir(n^) > 1 - 0(exp(-c(logL) 3 / 2 )) and finally 

\^ T (A)-ir(A)\ = 0(e~^ 3/2 ) 
for every event A e £l V) u, which implies W/J,^ — n\\ < l/(2eL 3 ) for L large enough. □ 



As a warm-up, we start by proving the easier Proposition 2. 

Proof of Proposition 2. Let £ e Sl^j/ satisfy (6.3). By monotonicity if rf & and £' e fi^t/ are 
such that £ < £' and 77 < 77', then 

P (at < L 10 : max (t) - 0°) > 2(log L) 3 / 2 ^ (6.7) 

< P' (at < L 10 : max (Vf'(t) - > 2 (log L) 3 ^ (6.8) 

where P' denotes the evolution with boundary condition 7/ (instead of 77). In particular, this 
is the case if we set rf x = r\ x + |_(3/2)(log L) 3 / 2 J and £' is sampled from the measure tt^(-\A), 
where A is the event A = A^^u = {(f) £ tt v ',u '■ 4>u > £,u}- From Theorem 1 (applied with 
e = a = 1/2) one sees that 

\\4(-\A)-4\\=o(e-< lo ^ 3/2 ) (6.9) 



for some c > 0. This is because 7/ is within distance ClogL from the plane ^ 3 / 2 )(\ogL) s / 2 ^ 

Definition 2) while £ is within distance (logL) 3 / 2 from the plane n n . Therefore, the probability 
in (6.7) is upper bounded by 

J -4{dg) P' (it < L w : max ($(t) - #?) > 2(logL) 3 / 2 ) + O (V^f^) 372 ) . (6.10) 

The initial condition £' in (6.10) is sampled from ir r ^, which is the invariant measure of the 

dynamics P', so that the distribution of cj£'(t) coincides with 7r^ at all later times. Via a union 
bound over times and recalling the relation between 77 and 77', the first term in (6.10) is upper 
bounded by 

4 (max (<j) x - 4%) > 2(logL) 3 / 2 ^ x 0(L 10 x L 2 ) (6.11) 
< vr^ (max {<j> x - 0£) > (l/2)(log L) 3/2 - lj X 0(L 10 x L 2 ) (6.12) 

which is of order exp(— c(logL) 3 / 2 ), see Theorem 1. The factor 0(L 10 x L 2 ) is just the average 
number of Markov chain moves within time L 10 , since there are order of I? lattice points in U. 
Similarly one bounds the probability that mm xe u(<t>x(t) — 4> x ) < — 2(logL) 3 / 2 for some t < L 10 
and claim (6.4) is proven. □ 
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Proof of Proposition 1. We prove only Eq. (6.1) since (6.2) is obtained essentially in the same 
way. Let W be a disk of radius 

p L = L x log L 

on the plane II" logL of slope n (cf. Definition 2, with C the same constant as in (2.4)) such that 
its projection V on the horizontal plane contains U and moreover the distance between dU and 
dV is at least pl/2 (recall that U has diameter L). 

Given u > 0, let C u be the spherical cap whose base is the disk W and whose height is u. The 
radius of curvature R is related to u and pi by 

(2R-u)u = p 2 L (6.13) 

and, since we will always work under the condition u <C pl *C -R, we have i? = p 2 / (2u)(l+o(l)). 
For a point u on the curved portion of the boundary of C u , let n„ be the normal at v directed 
towards the exterior of C u . It is clear that, if u <C pl, one has n„ = n + o(l); in particular, 
n v > with the convention of Definition 2. Finally the height (w.r.t. the horizontal plane) of 
the spherical cap at horizontal coordinates x G V is denoted by 

i/> u (x) = sup{z G M : (x (1) ,x (2) ,z) G C u }. 

We now define a sequence of spherical caps {C Un }^=o with constant base W, decreasing height 
u n and increasing radius of curvature R n . More precisely, let uq = 2L and M := 2L — (log L) 5 / 4 . 
Then we let u n = u n _i — 1 = 2L — n and (2R n — u n )u n = p\. For later purposes we also define 

*n = tn-l + #n (logL) 17/2 , t = 0. 

Remark 4. It is worth noting that R ~ L(logL) 2 /4, R M ~ L 2 /(2(log L) 3 / 4 ) and R n+ i/R n = 
l + o(l) uniformly in the whole range n = 0, . . . , M. 

Recalling that i? n = p L /(2u n )(l + o(l)), where o(l) is small uniformly in 1 < n < M, it is 
immediate to deduce that 

t M = pl(logL) 17 / 2 xo(^i]= 0(L 2 (logL) 23 / 2 ). (6.14) 

\n=l n / 

With this notation the key step is represented by the next Proposition. 

Proposition 4. There exists a positive constant d such that the following holds for L large enough. 
For every < n < M one has, with probability at least 1 — n exp(— c'(log -L) 3 / 2 ), 

4>x{t) - ^uS x )f or ever J x eU and every t G [t n , L 3 ]. (6.15) 

Since um <c (log L) 3 / 2 , Proposition 4 together with (6.14) imply the desired inequality (6.1). 

□ 

Proof of Proposition 4. We prove the claim by induction on n. For n = this is trivial since we 
chose no = 2L such as to guarantee that the maximal configuration A G £l v ,u is below the 
function U 3 x h-> i{j uo (x). 

Assume the claim for some n. For x G U define the event 

A x = {3t G [t n+1 ,L 3 ] : ^(t) > i> Un+ M} 

so we need to prove F(U x ^uA x ) < (n + 1) exp(— c/ (log L) 3 / 2 ). We have 

P(U X6 ^) < J2nA x ; A (s) < i, Un for every s G [t n , L 3 ]) + n e^'C ^) 372 (6.16) 

where we write <f> A (s) < ip Un to mean that <f> A (s) < il) Un (y) for every y G U. 
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Given x E U, consider the plane IT tangent to C Un at the point (x^, x^ 2 \ i^ Un ( x )) an d tne 
plane II' obtained by translating downwards II by (logL) 3//2 . The intersection of IT' with C Un 
is a disk W of radius 0(\/R~^ l (\ogL) 3 ^), whose projection on the horizontal plane we call Z. 
Let U C I? be such that U c Z and 5C7 is at distance of order 1 from dZ, so that of course 

diam{U) = 0{y/R^{\ogL) 3 ^). 

Let G O be the maximal monotone surface such that < ^ Un {x) for every x E U. Let 
P denote the law of the auxiliary monotone surface dynamics in U, starting at time t n from A^ 
and with boundary conditions A^~. By monotonicity and the definition of A^ we have 

W{A X ; <P h {s) < 4> Un for every s E [t n ,L 3 ]) < P (3i E [t n+1 ,L 3 ] such that ^ (t) > ^„ +1 (x)) . 
As in the proof of Theorem 2, Propositions 2 and 3 show that after time 

t n+1 - t n = R n+1 (log L) 17/2 > const x diam(J7) 2 (log L) 7 

the dynamics P has a variation distance of order exp(— c(log L) 3 / 2 ) for some c > from its 
equilibrium 7r^. ™ (recall that c_ log L < log R n < c + log L, cf. Remark 4) . Theorem 1 gives that 

Indeed, the point {x,ip Un+1 (x)) is at distance (1 + o(l))(log L) 3 / 2 from the plane IT containing 
the "planar" boundary condition A^. 

Putting everything together (plus a union bound over times t E [t n+ i,L 3 ] as in (6.11)) one 
gets 

F(A X ; c/) A {s) < ^ Un for every s E [t n , L 3 ]) = O (Vc(iogL) 3 / 2 ^ x o(L 3 x L 2 ) = O (e~% ( lo s L ) 3/2 ) . 
Finally, provided that we choose d = c/2, from (6.16) we get 

P(U. 6C /^)<(n+l)e- c '( lo s L ) 3/2 
(the union bound over x E U gives just an extra 0(L 2 )) which is the desired claim. □ 

7. SOS model: Proof of Theorem 3 

The proof of Theorem 3 is based on the following crucial results. The first is essentially an 
application of the results in [12] 2 , which can be formulated as follows. Let L and < h < L be 
fixed and define the set ri£ h of configurations £ £ ^L,h satisfying 

\Zi-m\ < v / Z(logL) K ; i = l,...,L, (7.1) 
where, as before, fji = [ih/(L + 1)J . 

Proposition 5. There exist a\ > and c > such that the following holds for every k > 1. For 
any L sufficiently large (how large depending on k), uniformly inO < h < L and £ E O^, ^: 

11/4 -tII <e- c(logL)2 , Vt > L 2 (logL) Ql . (7.2) 

Moreover, the same statement holds for the evolution constrained to stay above the wall f) (in this 
case of course one must require also that £ > fj). 



The results in [12] are stated for the discrete time chain, and thus a trivial overall factor of L must be taken into 
account when comparing our results with theirs. Moreover, [12] only deals with the case h = 0. 
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From the proof it will follow that one can actually choose a\ = 17 (with room to spare). The 
second crucial result deals with the relaxation of the extremal evolutions (the analogous result 
for monotone surfaces is Proposition 1). Let i] A ' + (t) denote the evolution of the maximal initial 
configuration A = L + h, with the wall constraint rji > fji, i = 1, . . . ,L. Similarly, let r/ v '~(i) 
denote the evolution of the minimal initial configuration V = — L, with the wall constraint 
m<Vi, i = l,--,L. 

Proposition 6. Let ot\ > be as in Proposition 5. For all L sufficiently large, and for some time 

7i = 0(1? (logL) 4+ai ): 

P(3iG{l,...,L}, 7/f' + (Ti) >f? i + v / I(logL) 5 ) <L~ 3 , (7.3) 

and similarly, 

P(3i€{l,...,L}, ^ V '"(Ti) < fji -JI (log L) 5 ) <L~ 3 . (7.4) 

Once Proposition 5 and Proposition 6 are established, it is an easy task to complete the proof 
of Theorem 3: 

Proof of Theorem 3. From a standard comparison estimate, cf. also (6.5) above, one has 

max ||/^ — vr 1 1 < CL 2 max (||/i A -vr||, ||^ -vr||) , 

where /i A , $ denote the law of evolutions r/ A (i), rj y (t) from maximal and minimal initial condi- 
tion respectively, with no wall constraint. Let G denote the event that |r/j — f)i\ < \< r L(logL) 5 for 
all i. Observe that for any event A, for t > T\ (where T\ is as in Proposition 6) 

\tf(A) - tt(A)\ < \F(r, A (t) G A ; ^(T,) G G) - „(A)\ + F^n) $ G) 
< max \\4_ Ti - tt\\ + 2P(t? a (T 1 ) $ G). 

By monotonicity one can couple the dynamics in such a way that r/ v '~(Ti) < ?? A (Ti) < ?? A,+ (Ti). 
Thus, from Proposition 6, P(7/ A (Ti) ^ G) < L~ 3 for all L large enough. On the other hand, 
from Proposition 5 one has max^g^ — M\ — e~ c ( logi ^ 2 , for T := L 2 (logL) Ql . Thus, taking 
t = T\ + T, it follows that ||/z A — vr 1 1 < 3L~ 3 , for all L large enough. The same bound, by 
symmetry, can be obtained for \\^ — vr 1 1 . Therefore, 

max \\uj — vr 1 1 < 3CL -1 , 
v 

which concludes the proof of Theorem 3, with any power a > a\ + 4. □ 

7.1. Proof of Proposition 6. This and the next subsection are devoted to the proof of Proposi- 
tion 6 (by symmetry it is sufficient to prove (7.3) only) assuming the validity of Proposition 5. 
The latter is proved in Section 7.3 below. We first recall a basic equilibrium estimate. 

Lemma 1. There exists some constant C > such that, uniformly in < h < L and H > y/Z log L: 

7r[3i = l,...,L: |& - fji\ > H] < Cexp {-C~ l mm{H 2 /L, H}). (7.5) 

Moreover, (7.5) continues to hold as it is for the interface conditioned to stay above the wall, as well 
as for the unbounded SOS model defined in (3.1). 

The proof of Lemma 1 can be obtained by standard large deviation arguments as in e.g. [12, 
Appendix C] . 

We give a detailed proof of Proposition 6 in the case h = only. The reader may verify that the 
same approach with minor modifications gives a proof for all < h < L. In order to underline 
the similarities between the proof of Proposition 6 with that of the companion result Proposition 
1 in the monotone surface case, we follow closely the notation introduced there. 
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Let pi = LlogL and, for u > 0, let C u C R 2 , denote the circular segment of height u and 
base a segment on the horizontal axis, with length 2pi and centered at L/2. We shall use 
the notation r/ £ C u for any configuration r/ such that < rn < t/j u (i),i = l,...,L where 
ipu(i) = sup{y G R : (i,y) £ C u }. Note that the radius of curvature i? of the circular segment C u 
satisfies (2R — u)u = p\. As in Section 6, in what follows we will always have u <C pl <■ R so 
that R = p 2 L /(2u)(l + o(l)). Define recursively no = 2L,to = and 

U n+1 =U n - (pl/u n ) l/3 (logL) 2 , t n+1 = t n + ( P l/u n f 3 (log L) 3+Ql 

and let M = min{n : u n < \fL(\og L) A }. Clearly, if R n denotes the radius of curvature of the 
circular segment C Un , then R n = p 2 L /(2u n ){\ + o(l)). 
Crucially, as in the monotone surface case, iju = 0(L 2 ). 

Lemma 2. For any L large enough tu < L 2 (logL) 4+Ql . 

Proof. We write 

M-l 

tu=J2 (piK) 4/3 (logL) 3+ ^ . (7.6) 

n=0 

Next observe that, by convexity, a 4 / 3 — 6 4 / 3 > 6 1 / 3 (a — 6) for all a > b > 0. Taking a = u n and 
6 = u n+ i, and using u n - u n+i = (p|/u n ) 1/3 (logL) 2 one has 

~ u% > uJUun - u n+1 ) > \L 2 l\\ogLfl* , 

where the last bound follows from u n+ \ > \u n which is easily seen to be implied by the assump- 
tion u n > \fL(\og L) 4 for all < n < M - 1. Therefore, 

M-2 

,.4/3 _ 4/3 V- / 4/3 _ 4/3 , 

m=n 

> u ^_ x + h M - n - l)L 2 / 3 (logL) 8 / 3 > i(M - n)L 2 / 3 (logL) 8 / 3 , 
the last bound following from u^^ 3 _ 1 > L 2 / 3 (logL) 8 / 3 . In conclusion, using this bound in (7.6), 

M-l 

t M < 2L 2 (logL) 3+Ql V < 2L 2 (logL) 4+ai , 

M — n 

n=0 

whenever L is large enough. □ 

Remark 5. As the careful reader has noticed, the length and time scales (u n , t n ) are quite different 
from their analogue in the monotone surface case. The main reason is the different order of mag- 
nitude of the maximal equilibrium fluctuations in the two models: log L versus \[L. However their 
value is determined by the common recipe which was described in Section 4. 

The key step in the proof of (7.3) is analogous to Proposition 4. 

Proposition 7. For any L large enough and for every < n < M, with probability at least 

l_ ne -(iogiO^ 

7] A (t) e C Un for every t £ [t n ,L 3 ]. 

If we now apply Proposition 7 with n = M and use the fact that = 0(L 2 (log L) 4+Ql ) and 
um < V~L(log L) 5 , we obtain the desired claim (7.3). □ 
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7.2. Proof of Proposition 7. As in the proof of Proposition 4 we proceed by induction in n < M. 
The initial step n = is obvious because the maximal configuration A is inside C UQ . Thus let us 
assume the statement true for n < M and let us prove it for n + 1. 

For 1 < i < L define the event 

A l = {3te [t n+1 ,L 3 ] : r£(t) > K+i(*)}- 
Using the inductive assumption we may write 

L 

F(uf =1 Ai) < J>(A; V A (s) G C Un for every s G [t n , L 3 ]) + iie"' logL ' 3/2 . (7.7) 

i=l 

Fix i = 1, . . . , L and consider the line L tangent to C Un at the point (i, tp Un (i)) and the line L' 
obtained by translating downwards (i.e. in the — y direction) L by 

2{u n - u n+1 ) = 2( / o|/n n ) 1 / 3 (logL) 2 . 

Let us denote by x_ , x + the horizontal coordinates of the leftmost and rightmost points of Uc\C Un 
and let / be the set of integers in [a;_,a; + ]. Clearly |/| = 0((p 2 L /u n ) 2 / 3 logL). 

Consider now the SOS dynamics in the interval / with boundary conditions equal to [Vv (*± )J 
at the left and right boundary i± of / and floor at zero height. This auxiliary evolution starts at 
time t n from the maximal configuration in the set of r/ G [0, L] 1 such that rn < ip Un (i) for 
any i E I. We denote by ¥' the law of this auxiliary chain. Observe that A^ n ) is within distance 
2(u n — u n+ \) = 0(^/]7[(log |/|) 3//2 ) from the line L' so that Proposition 5 will be applicable with 
k = 3/2. By monotonicity we have 

P (A; r] A (s) G C Un for every s G [t n , L 3 ]) < P' (it G [t n+1 , L 3 } such that ^ (t) > 4> Un+1 (t)) . 

Because of Proposition 5, after time |/| 2 (log < t n+ i — t n the dynamics P' has a variation 
distance of order exp(— c(logL) 2 ) from its equilibrium which we denote by n^. Since the 
distance between the point (i,ip Un+1 (i)) and the line L' is at least cy^7[(log |/|) 3,/2 for some 
c > 0, Lemma 1 gives that 

vrf } (m < W) > 1 - O (e-^I'D 3 ) > 1 - e -'0°^) 3 . 

As in the proof of Proposition 4, simple union bounds over i G [1, L] and t G [t n +i, L 3 ] give 

L 

^P(A i; v A (s) G C Un for every s G [t n ,L 3 ]) < L 5 e - c "( lo s L ) 2 < e -( lo s L ) 3/2 

i=l 

which finishes the proof of the inductive step. □ 

7.3. Proof of Proposition 5. We describe the main steps in the absence of the wall. At the end 
we will comment on the case with the wall. 

Let £ Ak (*)>£ Vk (*) denote the Glauber dynamics at time t starting from the maximal A K and 
minimal V K initial condition in Q£ h , respectively, and let //^ K , /j^ k denote their distribution. 
Notice that A K , V K are different from the maximal A and minimal V configuration in £li,h- F° r 
any £ G VL 1 ^ h , for any event Ac!!^ one has: 

4(A)-tt(A)= Y, <n)iA(A)-MA)). 
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Therefore, 



max ||/4 - vr || < Tr(Q LA \ Sl R Lh ) + max ||/4 - Mt II- 



By monotonicity, max^ >r)6 ^2 ; ||/4 — /^H < P(£ AK (i) / £ v,i (£))> where P(-) denotes the global 
monotone coupling (see Section 2.2). Using (7.5) to estimate h{VLl ,h \ ^1 h) f° r an y k > 1, we 
arrive at 



max wm 



tt|| < ce-^ L ^/ c + p(e AK (t) ^ e VK (*))■ 



Let now v + := 7r(- 1 £ > A K ) and := 7r(- 1 £ < V K ). Monotonicity implies that 

L L+h L L+h 

i=l j=~L i=l j=—L 

< 2L(L + h) (||^ + - ttII + - vrll) = 4L(L + h)\\tf - tt|| 

where /x^ is the law of the process at time t starting from the distribution v±. Above we used 
symmetry to conclude that \\^ + — ir\\ = \\fx^~ — ir\\. 

Since the event {£ > A K } is increasing, one has ir -< u + and the relative density v + (£)/-k(£) 
is an increasing function. The Peres-Winkler censoring argument [13] (see also Lemma 2.2 in 
[12]) proves that the same holds if we replace v + with . Thus the event F = (£) > 7r(£)} 
is increasing and \\^ + — ir\\ = ^i u t + {F) — tt(F). 

Next, we claim that, for some constants at\ < 17 and c > 0, for all t > L 2 (logL)" 1 and for 
any increasing event A, 

H u t + {A) -tt(A)< e-< los L ) 2 . (7.8) 

Clearly (7.8) finishes the proof of Proposition 5. In turn, the proof of (7.8) is adapted from the 
proof of a similar result in [12] (see Lemma 4.5 there). For the reader's convenience, we now 
describe its main steps. 

Step 1 . One first introduces an auxiliary parallel column dynamics which consists in replacing 
the single-site moves of the Glauber dynamics (3.1) with non-local moves as follows. At each in- 
teger time s = 1, 2, . . . firstly all even-numbered heights {^i}, i = 1,2,..., are re-sampled from 
their equilibrium distribution given the neighboring odd-numbered heights ^i+i, i = 0, 1, ... , 
and then viceversa with the role of even/odd columns reversed. By construction the column 
chain is ergodic and reversible w.r.t. it. Moreover, and that is the main reason for introducing 
it, Wilson's argument [18] applies to the new chain. Namely, by an almost exact computation, 
one can prove that under the column dynamics the maximal and minimal configuration in £lL,h 
couple in a time 0(L 2 logL) (see Section 3 in [12]). 

Step 2. Next, one relates the single site Glauber dynamics started from u + to the column dy- 
namics described above by means of the Peres-Winkler censoring idea. More precisely, given a 
time-lag T, one splits the time interval [0, t] into N = t/2T epochs each of length 2T (we assume 
for simplicity JVgN) and in each epoch one first censors (i.e. freezes) all the Poisson clocks at 
odd sites for the first half of the epoch and then all the even sites for the second half. If fi" + 
denotes the distribution of the censored chain, then the censoring inequality of [13] says that 
Vt + — ih + so that) f° r any increasing event A, 

tf(A) - ir(A) < ~tf(A) - ir (A) < ||^ + - tt||. 
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If the free parameter T is chosen in such a way that each epoch simulates very closely one step 
of the column chain then, thanks to Step 1, 

-A <e^ logL ) 2 

provided that N > L 2 (log L) 3 . 

As explained in [12], the overhead T introduced in the mixing time by the censoring is es- 
sentially the square of the maximum gradient |V£|oo := maxj — £j| that may arise in the 
censored dynamics. Since gradients of the equilibrium interface are i.i.d. random variables with 
exponential tail, conditioned to their sum being equal to h, well known estimates (see e.g. 
Lemma 1 above or [12, Appendix C]) imply that for any i = 0, . . . , L, uniformly in < h < L, 
and£> (logL) 2 : 

vr(|Ve|oo >£) <Ce- l ' c . 
Therefore, for any event A and starting from the distribution v := ir(-\A), the invariance of ir 
implies that 

fH (I V ^l- * l ) Z ^A) " ° W)' 

for I > (logL) 2 . A similar bound holds for the censored dynamics fi v t . The above observations 
lead to the following lemma 3 . 



Lemma 3 (see Lemma 4.2 in [12]). Let A be any increasing event, and consider the censored 
single-site dynamics started from v := ir(- | A). Let D := 1 + log(^) and choose T = D 2 {\ogLf 

and t = TL 2 (log L) 3 . Then, for some constant c > and for all L large enough 

II At - HI <e~ c ( logL ) 2 . 

Step 3. The above lemma cannot be applied directly to our case A = {£ > A K } because at 
time t = the maximum gradient is too large, of order \^L(logL) K , and therefore ir(A) = 
Q^ e -c\fL{\ogL) K ^ i s ver y sma ii -phis makes the corresponding overhead T too large. The solution 
to this problem proposed in [12] goes as follows (see also Lemma 4.5 and Corollary 4.6 there). 

Lemma 4. Fix k > 1. Let (logL) 2 < H < VL(logL) K and let A H = {£ i >f] i + H, i = l,... L}. 
Let uh '■= 7r(- 1 Ah) and let t = L 2 (log L) 16 . Then, for some c > 0, for L large enough, 

^(/)<^/ 2 (/) + e- c(logL)2 , (7.9) 
for any non-negative, increasing function f with \\f\\oo < L 

Proof. A sketch of the proof is as follows. Choose £ = H 2 /(logL) 2 and consider the Glauber 
dynamics on the enlarged state space £l£ t L,h of SOS interfaces {r]i}^^ e such that 

(rne[-L,L + h] ifie[l,L] 

\m > [ih/(L + 1)J if i G [-£, 0] U [L + 1, L + £] 

and with boundary condition 

V-(i + i) =-[(£+ l)h/(L + 1)J , VL+£+1 = [(L + £+ l)h/(L + l)j . 

Let 7r^ be the corresponding SOS reversible equilibrium measure and let [i t l,H be the distribu- 
tion at time t starting from u^h '■= 7rW(- 1 Ah)- Clearly the marginals of [i t l ' H and of ir^ on 



3 Actually in [12] the bound is given as o(l/L b ) for any b > but it is easily seen to hold in the form stated here. 
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£Il,h are stochastically larger than fj,^ H , and it respectively. Moreover, from standard local limit 
theorem estimates (see e.g. [12, Appendix C]) one has 

^){A H ) > J- e- CH2 / £ > e - c '^ 2 (7.10) 

for suitable constants C, C',7 > 0. Then, we apply Lemma 3 to the enlarged dynamics, with 
D = 0((log L) 2 ), and get that at time t = L 2 (log L) 6 , with b = 16, 

Thus, for any / as in the lemma, 

^ t H (f)<^' H (f)<^Hf) + e- cilosL)2 . 
Finally, let Eh = {?? G : {f/i < ff} n {t]l <h + H}}. Then, using the positivity of /, 

*® (/) < ttW (/ I ifc/a) + ttW (JS^/j) < tt(/ I ^ /2 ) + ttW (^ /2 ) . 

The Gaussian bounds of Lemma 1 show that 

^{E H/2 )<Ce- H2 / Ci <e- c( - l ^ L ^. 

In conclusion, the desired estimate (7.9) holds with t = L 2 (log L) 16 . □ 

Using the semigroup property, one can iterate the bound of Lemma 4 to go from the initial 
height Hq := \/L(k)gL) K to the final height H n := (logL) 2 , with n = O(logL) steps. Therefore, 
at time t = c'L 2 (log L) 17 , one has, for all / as in Lemma 4: 

^U)<^U) + e^^)\ (7.11) 

for some constant c', c" > 0. To finish the proof Proposition 5, cf. (7.8), observe that uh = v + 
and therefore, by (7.11), it is sufficient to establish 

th Hn (f) <^(/) + e- Cl ( logL )\ (7.12) 
for some c\ > 0, for all / as above. However, from the censoring inequality we know that 
lh Hn (/) — ih Hn (/)■ Moreover, a shift of the boundary height by 0((log V) 2 ) shows that 

^(^)>e- C(logi)2 . 

Therefore, (7.12) follows from Lemma 3, with D = 0((logL) 2 ), for any t > D 2 L 2 (log L) 12 = 
0(L 2 (logL) 16 ). This ends the proof of Proposition 5 without the wall, with the choice of the 
constant a\ = 17. 

When the wall is present, the symmetry between the maximal and minimal configuration 
breaks down and one has the new problem of proving convergence to equilibrium within the 
correct time scale starting from the minimal configuration fj. This problem has been solved in 
[12] when h = but the same proof applies to h € [0, L\. □ 

8. SOS model: Proof of Theorem 4 

Let gap(L, H) denote the spectral gap of the SOS model with zero boundary conditions, floor 
at —H and ceiling at +H. Similarly, let gap tu (L,i7) denote the same quantity for the system 
with the wall, i.e. the spectral gap of the SOS model with zero boundary conditions, floor at 
and ceiling at +H. 

Proposition 8. There exists c > such that for all L e N sufficiently large, and for all H e N: 

gap(L,#) >cL- 2 (logL)- Q , (8.1) 
where a > is the same constant appearing in Theorem 3. The same bound holds for gap"'(L, H). 
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The first observation is that Proposition 8 immediately implies Theorem 4. Indeed, this fol- 
lows from an elementary approximation since the bound (8.1) is uniform in H: both £ (/, /) and 
Var(/) can be seen as limit as H — > oo of the corresponding quantities for the model with floor 
and ceiling at ±H, so that for fixed L one has gap(L) > liminf^^oo gap(L, H). 

8.1. Proof of Proposition 8. To prove Proposition 8 we use a recursive approach together with 
a powerful idea borrowed from the mathematical theory of kinetically constrained spin systems 
(see section 4 in [2] and the proof of Proposition 9 below). We give the details of the proof of 
the lower bound on gap w (L, H), and will later illustrate the minor modifications necessary to 
cover the case without the wall; see Section 8.3 below. 

For any rectangle A in I?, with base £ A (horizontal side) and height h\ (vertical side), we 
write 7ta for the equilibrium measure of the SOS model in the rectangle A with zero boundary 
conditions at the endpoints of the bottom horizontal side of length £\ (the lower horizontal side 
of the rectangle has vertical coordinate zero) . We call Q\ the set of all allowed configurations, 
i.e. £1 A = {0, . . . , /ja} £a . Let also 

(8 ' 2) 

where we use the notation £ A , Yai\ for the Dirichlet form and the variance associated to ir A , 
and / ranges over all functions such that VarA(/) / 0. Note that, with this notation, one has 

gap(A) = gap w (£ A ,h A ). Define 

7(L,n)= max gap(A) -1 , 

where F„ l is the class of all rectangles A in Z 2 such that £\ < L and h\ < (3/2) n . The crucial 
recursive estimate reads as follows. 

Proposition 9. Let n = n (L) £ N be such that (3/2) no < L < (3/2) no+1 . Then for any L large 
enough, for any n > uq, 

7(L,n)<(l + <S n ) 7 (L J n-l) (8.3) 

where 8 n = (3/2)~ n / 6 . 

Let us finish the proof of Proposition 8, assuming the validity of Proposition 9. Iterating the 
bound (8.3) one has 

j(L,n) <Cj(L,n ), (8.4) 

and 



C = ] [ (1 + S n ) < oo, 



n=no 



uniformly in L. Proposition 8 now follows from (8.4) and the following lemma. 
Lemma 5. There exists C < oo such that for any L large enough 

7 (L,n ) <CL 2 (log L) a , 

with no = no(L) as in Proposition 9. 

Proof. Note that, thanks to the straightforward bound gap -1 < T m i x (which holds for any re- 
versible Markov chain [10]), if h A < (3/2)™° < C£ A < CL for some C > 0, from Theorem 3 one 
has 

gap(A)^ 1 = 0(L 2 (logLn. 
(Strictly speaking Theorem 3 is stated only for h\ < £ A , but the proof works as it is as soon 
as h A < Cl A for some C). Therefore, the same bound applies to every A e ^ no ,L such that 
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£a > L/8. Thus, one has to estimate only max{gap(A) _1 , A e ^ no ,L '■ £\ < L/8}. In order to 
deal with this term, define f(L) = j(L, no), for all L and no = no(L) as above. From Proposition 
9, one has 

max{gap(A)- 1 , A G F no , L : i k < L/8} < d /(L/8) , 
where the constant C\ = IIn=no-3^ + $n) < 2 for no large enough. In conclusion, 

f(L) < max{2/(L/8),CL 2 (logL) a } < 2f(L/8) + CL 2 (log L) a . 

A simple induction now proves that f(L) < C'L 2 (log L) a for some new constant C. □ 

8.2. Proof of Proposition 9. Take n e N and A £ F n> £. Without loss of generality we write 
A = [1,1a] x [0, h\]. We need to show that 

gap(A)- 1 < (1 + 5 n )j(L,n- 1). 

We use a geometric construction analogous to that of [1, Proposition 3.2]. If h k < (3/2) n ~ 1 
then actually A e F„_i ;L and therefore gap(A)^ 1 < -y(L,n - 1). If instead h\ > (3/2)™" 1 , then 
it is easily seen that we can write A = Ai U A2 where Ai, A2 are such that: 

(a) ^ e F n _i jL , i = 1,2; 

(b) Aj, i = 1,2, has the same base [1,1a] of A; 

(c) the overlap rectangle / = Aifl A 2 / has base length £\ and height at least A n := (3/2) 3n / 4 . 
Moreover, there are at least s n := (3/2) n / 5 such decompositions {(A^, A^)} s ™ with the prop- 
erty that the overlap rectangles = A^ n Ag are disjoint, i.e. 1^ n 1^ = 0, for all i 7^ j. 

Next, fix one of the s n decompositions mentioned above, say Ai = [1,£a] x [0,/iaJ and 
A2 = [1, £a] x [h-Ai — \I\,hA], where h-A = h,A 1 + ^a 2 — hi, with /i/ > A n denoting the height of 
the overlap rectangle /. Consider the distribution 7r Aa corresponding to the equilibrium measure 
in the region A2 conditioned to the value of the configuration r/ in the region A \ A2 = [1,-^a] x 
[0, ZiAi — hi]. Let 

C(r,) := {i £ [1,£ A ] : rji > h Al - hi}. 

Note that, if C(rj) 7^ 0, then C(rf) is the disjoint union of intervals B\,B2,..., B^, and one has 
that 7r A is a product measure ®f =1 7Ti^ on the rectangles ii.; = Bi x [/i Al — /tj + 1, Iia], such that 
each 7rj^ is the SOS equilibrium measure with floor at height hA 1 — hi + 1 and ceiling at height 
/ia- If, on the other hand, C(n) = 0, then 7/ is entirely contained inside A \ A2, and 7r Aa is trivial, 
i.e. it gives full mass to the empty configuration in the region A2. 

Define a constrained block dynamics as follows. With rate one the current configuration 77 
is re-sampled inside A2 from the distribution tt'^ described above. Moreover, if n c Ai, that is 
Vi < ^Ai for all i = 1, ... , £a, then with rate one, the current configuration 77 is re-sampled inside 
Ai from the distribution 7r Al (defined just before (8.2)). This process is kinetically constrained 
in that the region Ai is updated only if n C Ai. 

It is not hard to check that the above dynamics is reversible w.r.t. the SOS equilibrium measure 
7ta, and that its Dirichlet form is equal to: 

£block(/,/) = vr A (l{^cAi}Var Al (/)+Var A2 (/)^ , 

where Var Al stands for the variance w.r.t. 7ta 13 while Var Aa stands for the variance w.r.t. 7r\ 2 . 
Notice that, even starting from the maximal configuration A in A, if we first update A2 and then 
Ai, with very high probability we reach the distribution tta 1 ■ Indeed, thanks to the fact that hi > 
A n » \J~£a, the first update in A2 will produce with very high probability a new configuration 
entirely contained in Ai. In turn, since \\tta 1 —tta\\ is very small (because ^i Al » y/iX), the above 
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example suggests that the spectral gap of the constrained block dynamics should be very close 
to one. This is quantified in the next lemma. 

Lemma 6. For all L large enough, and n > no(L), for all functions f: 

Var A (/)< (l + e-( 3 / 2 )" /4 )f block (/,/). 

Assuming the validity of Lemma 6, the proof of Proposition 9 can be completed as follows. 
From (3.6) and (8.2), for all / : n A R: 

(^cA^VarAiC/ - )) 

1 lA 

< 7 (I, n- l)-^^ A [ K , + (V ii+ /) 2 l te < /lAi} + K ,_(V i ,_/) 2 l te < /lAi} ] . (8.5) 
i=l 

Similarly, 

7r A (Var A2 (/)) 

} ] • (8.6) 

z i=i 

From Lemma 6, (8.5) and (8.6) we have 

Var A (/) < (1 + e-^ n/4 ) 7 (L, n - 1) (S A (/, /) + /)) (8-7) 
where £\(f,f) is the usual Dirichlet form in A, while 

1 ' A 

S {(f>f) '■= 2 5Z 7rA [ Pi '+( Vi > + ^) 2l { /l A 1 -^<%<^A 1 } +^-( V * -/) 2l {^A 1 -^/<%<^A 1 }] • 
1=1 

Since (8.7) is valid for every one of the s n choices of the rectangles (A^ , A2 )*™ 13 one can 
average this estimate to obtain, with = A^ n Ag : 

Var A (/) < (1 + e -W n/A ) 7 (L, n - 1) (£ A (f, f) + s" 1 £ (/, /)) . 

Since the overlaps 1^ are disjoint, one has the obvious bound Ylj=i ^A^ (/> /) — £a(/>/)- It 
follows that 

7 (L, n) < (1 + e -( 3 / 2 )" /4 ) (l + S ^) 7 (L, n - 1). 

By construction, s" 1 = 0((3/2)" n / 5 ), so that (1 + e^( 3 / 2 )" /4 )(l + s" 1 ) < 1 + 5 n for all n large 
enough, with <5„ = (3/2) _n / 6 . This concludes the proof of Proposition 9. 

Proof of Lemma 6. The proof is similar to that of [2, Proposition 4.4]. Let G = {77 G S7 A : 77 C 
Ai}. Note that the event G is very likely to occur under the equilibrium 7r A . In particular, from 
(7.5), one has 

e ■= tt a (G c ) < C exp (-C _:l (3/2) n ). (8.8) 

Moreover, G is very likely to occur even under the equilibrium 7r Aa , uniformly in 77 e J7 A . Indeed, 
by monotonicity the smallest value of 7r Aa (G) is achieved at configurations 77 G f2 A such that 
?7i > ^Ai — /i/ for all i = 1, . . . ,£\ and in that case tt\ (G) is the 7r A / -probability in the rectangle 
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A' 2 := [1, £\] x [0, h\ 2 ] that the height does not exceed hi at any point. Therefore, using hj/£^ > 
A 2 (3/2)" n = (3/2) n / 2 and the bounds (7.5) one has 

5 := max7r A (G c ) < C exp (-C" 1 (3/2) n / 2 ). (8.9) 



On the other hand, reasoning as in (7.10) (and using h\ 2 » 5 also satisfies 



and therefore in particular 6 > y/e/{\ — e) for all n large enough. 

The infinitesimal generator of the block dynamics acts on functions / : Oa h> E as 

Adock/fa) = IgW (tt Ai (/) - /fa)) + <(/) - /(r?) . 

Let A be the spectral gap of the block dynamics and let / be an eigenfunction of £t>iock with 
eigenvalue — A. Thus, / must satisfy, for all rj E Oa- 

(1 - \)f(rj) = 1 G ( V ) (vr Al (/) - f( V )) + <(/). (8.10) 

If we average both sides w.r.t. tt\ l = vta(- | G) we get 

(1 - A)tt Ai (/) = tt Ai «(/)) = vr A « 2 (/) | G) = (1 - 6)" 1 Cov A (l G , <(/)) (8.11) 

where Cova(- , •) denotes the covariance w.r.t. 7r A , and we have used 7ta(7t a (/)) = 7ta(/) = 0. 

Assume now 1 — A > 0, otherwise there is nothing to be proved. Schwarz' inequality gives 
Cov A (l G , ttX 2 (/)) < Var A (l G ) 1/2 |KA 2 (/)||oo. Therefore, from (8.11), one has 



K Al (/)| < TfrW ^ IKa 2 (/)IU < T^^y IKa 2 (/)IU, (8.12) 



where we use 8 > \fe]Jl — e), for n large enough and ||7r A (/)||oo = su Pr; l 7r A 2 (-^)l- ® n tne otner 
hand we can rewrite (8.10) as 



1-A + 1 g (t?) 1-A + 1 g (t7) 
Therefore, applying 7r A2 to both sides, we get 

^ = 2-A + 2-A + i^A • 

Then, using (8.12), one has 

IKa 2 (/)IIoo < (j^ C^Tx + + 2^a) (/)l|o ° =: X(<5 ' A) l^miloc. 

The latter bound is possible only if: (a) 7r A (/) = or (b) x(<5, A) > 1. It is easy to exclude (a). 
Suppose in fact that 7r A (/) = 0. Then also 7ta 1 (/) = by (8.12), and (8.10) would reduce to 
(1 — A + 1 G )/ = which is possible iff / = because of the assumption 1 — A > 0. On the other 
hand it is straightforward to check that x(<5, A) > 1 implies that 1 — A = 0(V6) < e~( 3 / 2 ) n 4 for 
all n large enough. The proof of the Lemma is complete. □ 
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8.3. Spectral gap lower bound without the wall. The proof of Proposition 8 in the case with- 
out the wall is very similar. The analogue of Proposition 9 is described as follows. Given a 
rectangle A = [1,£a] x [— h\,h\], write gap(A) = gap(-^A,^A) for the spectral gap of the SOS 
model in A with zero boundary conditions at and £\ + 1. Set 

P(L,ri) = gap(A) - " 1 , 

where G n> i stands for the set of rectangles A = [1,£a] x [— ^A, ^a] such that l\ < L and h\ < 
(3/2)". The decompositions presented in items (a)-(b)-(c) in Section 8.2 can now be obtained 
as follows. Let Ai = [1,£\] x [—hA 1 ,h\ 1 ] and A2 be given by the union of two rectangles 
A2 = Ag U A\, where the bottom rectangle is h\ = [1,£a] x [—h\,—h\ + h\ 2 ] while the top 
rectangle is A 2 = [1,£\] x [Ha — ^a 2 j^a]- Note that if h\ 1 + Ha 2 > ^a then there are two 
overlap regions I b ,P. As in Section 8.2, if A e G n ,L is such that Jia > (3/2) n_1 then one 
can find at least s n decompositions of the form A = Ai U A2 where Aj, i = 1, 2 are as above 
with /i Al + Ha 2 - Ha > 2A n , with A x e G n _ ljL and A^,A^ e F n _ 1L , and such that all the 
overlaps I = P U I b corresponding to distinct decompositions are disjoint. Let tta, tta 1 be now 
the equilibrium distribution in A, Ai with zero boundary conditions and ir\ be the distribution 
in the region A2 conditioned to the value of rj G A \ A2. With this notation, it is not hard to 
check that the result of Lemma 6 remains true as it stands. Then, the same argument of Section 
8.2 proves that 

p(L,n) < (l + <y n )fen-l), 
where we define (3(L,n) := max{/3(L, n), j(L, n)}. Since we know the results for the con- 
stants j(L, n) (cf. Proposition 9), it is now simple to infer the desired conclusion: f3(L, n) < 
CL 2 (log L) a . This is sufficient to end the proof in the case without the wall. 
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